*     discharge influence function at CZ for areasink
*     with corners CZ1 CZ2 CZ3 CZ4i in counterclockwise order
      COMPLEX*16 FUNCTION CFSAQ0(CZ, CZ1, CZ2, CZ3, CZ4)
      IMPLICIT NONE
      COMPLEX*16 CZ, CZ1, CZ2, CZ3, CZ4

      COMPLEX*16       CZ2M1I, czt3, czt4, cal1, cal2, cal3,
     &                 CZT, CZTSQ, C1ZT, CZTM1, C1ZTM1, CZTMZ3,
     &                 C1ZMZ3, CZMZ3S, CZTMZ4, C1ZMZ4,
     &                 CF1, CF2, CF3, CF4, CFLG1, CFLG2, CFLG3, CFLG4,
     &                 CFLGT, CFD1, CFD2, CFD3, CFD4, CFD5, CFD, CFSAW
      DOUBLE PRECISION RKSI3, RKSI4, RET3, RET4, ra1, ra2, RL, rcof,
     &                 RKSI, RTKSI3, RKSI3S, RKSIM3

      DOUBLE PRECISION RO2PI
      COMPLEX*16       CI
      PARAMETER (RO2PI=.159154943,CI=(.0,1.))


      CZ2M1I=1./(CZ2-CZ1)
      czt3=(CZ3-CZ1)*CZ2M1I
      czt4=(CZ4-CZ1)*CZ2M1I
      RKSI3=REAL(CZT3)
      RKSI4=REAL(CZT4)
      RET3=AIMAG(CZT3)
      RET4=AIMAG(CZT4)
      cal1=(RKSI3-1.)/(CZT3-1.)
      cal2=(RKSI3-RKSI4)/(CZT3-CZT4)
      cal3=RKSI4/CZT4
      ra1=(1.-RKSI3)*RET3
      ra2=(RKSI3-RKSI4)*(RET3+RET4)+RKSI4*RET4
      RL=ABS(CZ2-CZ1)
      rcof=.5*RL*RL*RO2PI

      CZT=(CZ-cz1)*cz2m1i
      CZTSQ=CZT*CZT
      RKSI=REAL(CZT)
      C1ZT=1./CZT
      CZTM1=CZT-1.
      C1ZTM1=1./CZTM1
      CZTMZ3=CZT-czt3
      C1ZMZ3=1./CZTMZ3
      CZMZ3S=CZTMZ3*CZTMZ3
      CZTMZ4=CZT-czt4
      C1ZMZ4=1./CZTMZ4
      RTKSI3=REAL(czt3)
      RKSI3S=RTKSI3*RTKSI3
      RKSIM3=RKSI-RTKSI3
      CF1=2.*CI*(-CZT+RKSI)
      CF2=-2.*CI*(cal1*CZTMZ3-RKSIM3)
      CF3=2.*CI*(cal3*CZT-RKSI)
      CF4=-2.*CI*(cal2*CZTMZ3-RKSIM3)
      CFLG1=CF1*LOG(CZTM1*C1ZT)
      CFLG2=CF2*LOG(CZTMZ3*C1ZTM1)
      CFLG3=CF3*LOG(CZTMZ4*C1ZT)
      CFLG4=CF4*LOG(CZTMZ4*C1ZMZ3)
      CFLGT=CFLG1+CFLG2+CFLG3+CFLG4
      CFD1=-(CI*CZTSQ+ra1)*C1ZT*C1ZTM1
      CFD2=-CI*(cal1*CZMZ3S+2*CZT*RTKSI3-RKSI3S)
      CFD2=CFD2*(czt3-1.)*C1ZTM1*C1ZMZ3
      CFD3=-(-CI*cal3*CZTSQ+ra2)*czt4
      CFD3=CFD3*C1ZT*C1ZMZ4
      CFD4=-CI*(cal2*CZMZ3S+2.*CZT*RTKSI3-RKSI3S)
      CFD4=CFD4*(czt4-czt3)*C1ZMZ3*C1ZMZ4
      CFD5=-(ra1+ra2)*C1ZT
      CFD=CFD1+CFD2+CFD3+CFD4+CFD5
      CFSAW=(CFLGT+CFD)*rcof*cz2m1i
      CFSAQ0=CONJG(CFSAW)

      RETURN
      END
